Astronomy &: Astrophysics manuscript no. ruplim 


© ESO 2012 


October 2, 2012 





Constraints on off-axis jets from stellar tidal disruption flares 

Sjoert van Velzen 1 , Dale A. Frail 2 , Elmar Kording , and Heino Falcke 1 ' 3 ' 4 

1 IMAPP, Radboud University, P.O. Box 9010, 6500 GL Nijmegen, The Netherlands 
e-mail: s.vanvelzen@astro.ru.nl 

2 National Radio Astronomy Observatory, Socorro, NM, USA 

3 ASTRON, Dwingeloo, The Netherlands 

4 Max-Planck-Institut fur Radioastronomie Bonn, Germany 

Preprint online version: October 2, 2012 

Abstract 

Context. Many decades of observations of active galactic nuclei (AGN) and X-ray binaries have shown that relativistic 
jets are ubiquitous when compact objects accrete. One could therefore anticipate the launch of a jet after a star is 
disrupted and accreted by a massive black hole. This birth of a relativistic jet may have been observed recently in two 
stellar tidal disruption flares (TDFs), which were discovered in gamma-rays by Swift. Yet no transient radio emission 
has been detected from the tens of TDF candidates that were discovered at optical to soft X-ray frequencies. Because 
the sample that was followed-up at radio frequencies is small, the non-detections can be explained by Doppler boosting, 
which reduces the jet flux for off-axis observers. And since the existing followup observation are mostly within ~ 10 
months of the discovery, the non-detections can also be due to a delay of the radio emission with respect to the time 
of disruption. 

Aims. We wish to test the conjecture that all TDFs launch jets. 

Methods. We present 5 GHz follow-up observations with the Jansky VLA of seven known TDFs, doubling the number 
of radio observations of these events. To avoid missing delayed jet emission, our observations probe 1-8 years since the 
estimated time of disruption. 

Results. None of the sources are detected, with very deep upper limits at the 10 micro Jansky level. These observations 
rule out the hypothesis that these TDFs launched jets similar to radio-loud quasars. We also constrain the possibility 
that the flares hosted a jet identical to Sw 1644+57, the first and best-sampled relativistic TDF. 

Conclusions. We thus obtain evidence for a dichotomy in the stellar tidal disruption population, implying that the jet 
launching mechanism is sensitive to the parameters of the disruption. 



1. Introduction 

The disruption of a star by a massive black hole leads to 
arguably the most spectacular form of accretion onto these 
compact objects. The stellar debris that remains bound af- 
ter the disruption returns to the black hole at a rate that 
initially can exceed the Eddington limit (MEdd) by many 
orders of magnitude. While the rapid i~ 5 / 3 decline of this 
fallback rate (Rees 1988) implies that 1% of MEdd is typi- 
cally reached within a few to ten years. A tidal disruption 
flare (TDF) may thus be used to sample different modes 
of accretion (e.g., Abramowicz & Fragile 2011) for a sin- 
gle supermassive black hole. Considerable effort is needed 
to simulate the dynamics of the disruption (e.g., Nolthcnius 
& Katz 1982; Evans & Kochanek 1989; Rosswog et al. 2009; 
Guillochon & Ramirez-Ruiz 2012) and to estimate the re- 
sulting optical to X-ray light curve of the flare (e.g., Loeb 
& Ulmer 1997; Bogdanovic et al. 2004; Strubbe & Quataert 
2009; Lodato & Rossi 2011). Efficient detection to obtain 
large samples TDFs is much anticipated, as this will allow, 
for a example, a study of the demographics of dormant 
black holes beyond the local universe (Frank & Rees 1976; 
Lidskii & Ozernoi 1979). 

Tens of (candidate) stellar tidal disruption events have 
been found by searching for flares in soft X-ray (Komossa 
& Bade 1999; Grupe et al. 1999; Komossa & Greiner 1999; 



Greiner et al. 2000; Esquej et al. 2008; Maksym et al. 2010; 
Lin et al. 2011; Saxton et al. 2012), UV (Gezari et al. 2006, 
2008, 2009, 2012), or optical surveys (van Velzen et al. 2011; 
Drake et al. 2011; Cenko et al. 2012a), or based on spec- 
tra with extreme coronal lines (Komossa et al. 2008; Wang 
et al. 2012). None of these thermal flares are associated with 
a radio transient, but only a handful have been followed up 
at this frequency. The only tidal disruption candidates with 
a detected transient radio counterpart are those discovered 
in 7-rays by Swift: Sw 1644+57 (Bloom et al. 2011; Burrows 
et al. 2011; Levan et al. 2011; Zauderer et al. 2011) and 
Sw 2058+05 (Cenko et al. 2012b). Since the radio and X- 
ray emission of these two events most likely originates from 
a relativistic jetted outflow, they are often referred to as 
relativistic TDFs. In this paper we shall refer to the other 
class of TDFs as 'thermal', since they are all discovered at 
optical to soft X-ray frequencies. 

One is left to wonder why the two TDFs discovered 
with Swift are the only events with evidence for a newly- 
born jet. Interpreting this as a radio-loud/radio-quiet di- 
chotomy, similar to the devision of radio- loudness in quasars 
(Kcllermann et al. 1989; Falcke et al. 1996; Sikora et al. 
2007), would require that the tidal disruption jet launch- 
ing mechanism is sensitive to the properties of the disrup- 
tion (e.g., mass ratio, impact parameter, orientation of the 
orbit of the star with respect to the black hole spin, or 
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circumnuclcar environment). The explanation that quasars 
spend only a fraction of their time as radio-loud objects, 
similar to the jets in the 'hard intermediate state' of X- 
ray binaries (e.g., Kording et al. 2006), does not apply to 
tidal disruptions because their accretion rate is not con- 
stant. On the other hand, based on the observation of a 
fundamental plane of black hole accretion, a universal scal- 
ing law for the non-thermal emission across the entire black 
hole mass range (Merloni, Heinz, & di Matteo 2003; Falcke, 
Kording, & Markoff 2004), and the abundance of jets low 
luminosity AGN (Nagar et al. 2000) and X-ray binaries or 
microquasars (Mirabel & Rodriguez 1999; Fender 2001), 
one may postulate that all stellar tidal disruptions launch 
jets. Likewise, Miller & Giiltekin (2011) argue that the fun- 
damental plane can be used to estimate the black hole mass 
of a TDF (if X-ray and radio observations of the flare are 
available) . 

If all stellar disruptions are indeed companied by rela- 
tivistic outflow, the current upper limits on the radio flux 
of the thermal TDFs can be explained by the orientation of 
this jet, which can dramatically reduce the flux due to rela- 
tivistic Doppler boosting, and by a delay of the radio emis- 
sion of the jet with respect the time of disruption. However, 
the number of TDFs that have been followed-up at radio 
frequencies is currently not sufficient to test this unification 
based on viewing angle. 

Recent advances in the hardware of the Very Large 
Array (VLA) have made it possible to obtain very deep 
radio observations of stellar tidal disruptions in a relatively 
short amount of time. To use this opportunity, we selected 
all thermal stellar tidal disruptions that occurred after 2004 
for follow-up observations. These observations more than 
double the number of TDFs with deep radio observations. 
And because our radio observations span a wide range of 
times since the disruption, we can, for the first time, test 
the hypothesis that all stellar tidal disruptions launch jets. 

The rest of paper is organized as follows. In sec. 2 we 
present two different tidal disruption (TD) jet models and 
compute off-axis light curves. In sec. 3 we discuss the radio 
observations and sample selection. We use these observa- 
tions to constrain the jet models in sec. 4 and we close with 
a discussion in sec. 5. 



2. Tidal disruption jet models 

To be able to interpret our radio observations, we need a 
model that describes how the radio emission in jets of ac- 
creting objects is typically made. In this section we there- 
fore review two models of tidal disruption jets and we 
present off-axis light curves for these models. We have di- 
vided the models into two classes 1 based on the origin of 
the emitting particles: external or internal. In both models, 
some fraction of the accretion power ends up in the jet and 
the emission mechanism is synchrotron radiation. 



2.1. External model: off-axis light curves for Sw 1644+57 

The external model of radio emission from TD jets was 
first presented by Giannios & Metzger (2011) and further 
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Figure 1. The observed light curve of Sw 1644+57 (open 
symbols), with the predicted late-time light curve (dashed 
lines) for a total jet energy of Ej = 10 52 erg (Berger et al. 
2012). We show the estimated 5 GHz light curve of different 
off-axis observers, Eq. 2, assuming that the Lorentz factor 
of the jet decreases with Tj oc t~ - 2 , as inferred by Berger 
et al. (2012). We modified the extrapolated light curve to 
match a Sedov- Taylor solution, Lj oc i -1 ' 2 , when Tj < 2. 
The 2-cr upper limits on the radio flux of seven other TDF 
candidates (Table 2) are shown with black triangles (we 
scaled these limits to the redshift of Sw 1644+57, see sec. 
4.1). 



developed by Metzger, Giannios, & Mimica (2012). Shock 
interaction between the jet and the gaseous circumnuclear 
medium powers the emission, similar to afterglow models 
of gamma-ray bursts (e.g., Sari & Piran 1995). 

The external model has been applied to the radio light 
curve of the relativistic TDF Sw 1644+57 (Metzger et al. 
2012; Berger et al. 2012), we show the fit and predicted late- 
time light curve in Fig. 1. We note that this fit requires an 
continuous increase of the isotropic jet power during the 
first year of observations. 

The scaling of the synchrotron peak and self absorption 
frequency in the Metzger et al. (2012) model of Sw 1644+57 
are based on spherical expansion of an ultra-relativistic 
shell and thus require OjTj < 1 (with 6j, Tj the jet opening 
angle and Lorentz factor, respectively), plus an on-axis ob- 
server i\ < l/Tj(t = 0); both requirements are supported 
by the observed radio light curve (Metzger et al. 2012). 

To compute the light curve for an off-axis observer, we 
first boost the observed on-axis flux Fi(v) into the jet rest- 
frame 



L j {v) = di8 s l - a F 1 {v) 



(1) 



1 Other models of TD jets (Lei & Zhang 2011; Krolik & Piran 
2011; De Colle et al. 2012), are not discussed here since these 
make no predictions for off-axis light curves at a given observer 
frequency. 



(e.g., Lind & Blandford 1985; Jester 2008) here we intro- 
duced the Doppler factor for the on-axis observer S\ = 
\Tj(l — /3j cos?i)] _1 with (3j = Vj/c, a is the spectral index 
defined as F(v) cx v a , and is the luminosity distance. 
Next, we transform the jet luminosity to the off-axis ob- 
server using a different Doppler factor, 82 . If the size of the 
emitting region is small compared to the distance to the 
black hole, the time delay due to the geometrical separa- 
tion of the synchrotron peak with frequency can be ignored, 
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and we can estimate the flux for an observer sitting at i^: 



D3-13 (5 GHz) 



St 



l-ft(t) 

1 — (3j (t) COS 12 



3-a(t) 



F x {t,v). (2) 



Here t is measured in the observer-frame, a(t) is obtained 
from the light curve, and we used cosii ~ 1 — TJ 2 /2 ~ 
1 to simplify the equation. If Sw 1644+57 was indeed a 
rclativistic outflow that we observed on-axis, the light curve 
for off-axis observers depends only Tj(t) and i 2 . The latter 
is a free-parameter, which we shall constrain by our follow- 
up observations (in sec. 4). 

The external model of TD jets (Metzger et al. 2012), 
applied to Sw 1644+57 by Berger et al. (2012), is used to 
obtain Tj(t) and the light curve beyond the last published 
radio observation of Sw 1644+57 (about one year after the 
Swift trigger). The off-axis light curve that is derived here 
may thus be viewed as a test for this model. We consider 
two scenarios. First, we set T(t > 1 yr) = 2, and use the 
extrapolated light curve presented in Berger et al. (2012). 
This constant Lorentz factor is required because the extrap- 
olated light curve is no longer valid when the jet slows down 
to mildly relativistic speed and lateral expansion becomes 
important. This happens at Tj < 2 (Zhang & MacFadyen 
2009). We also consider a decreasing jet velocity Tj cx i~°- 2 
(as inferred for Sw 1644+57), but modify the extrapolated 
light curve to match the non-relativistic Sedov- Taylor evo- 
lution Lj cx t~ 12 (e.g., Granot et al. 1999; Leventis et al. 
2012) when Tj < 2. We show the light curves in Fig. 1. 



2.2. Internal model: off-axis light curves for known TDFs 

The internal model of radio emission from TD jets was 
first presented in van Velzen, Falcke, & Farrar (2010) and 
further developed by van Velzen, Kording, & Falcke (2011). 
The model is based on the simple idea of jet-disk coupling 
(Rawlings & Saunders 1991; Falcke & Biermann 1995): a 
constant fraction of the accretion luminosity (Lj) is fed into 
the jet, Qj — qjLd- The conversion from jet power (Qj) 
to radio luminosity (Lj) follows by assuming equipartition 
between the energy in relativistic particles and magnetic 
fields, and has been calibrated using observations of AGN 
(Falcke et al. 1995; Willott et al. 1999; Kording et al. 2008). 

Stellar mass black holes show rapid switches from radio- 
loud (qj = 0.2) to radio-quiet (qj < 0.002) coupling as the 
accretion rate increases from sub-Eddington to (near) the 
Eddington limit (Fender, Belloni, & Gallo 2004). Motivated 
by the growing evidence that accretion onto super-massive 
black holes can also divided into these two modes (e.g., Ho 
1999; Ghisellini & Celotti 2001; Falcke et al. 2004; Kording 
et al. 2006; Best & Heckman 2012; Plotkin et al. 2012), 
we considered the following three scenarios for the jet-disk 
coupling in tidal disruptions: 




all times (a) 
M(t) > 2%M Edd (b) 

t < ifallback (c) 



(3) 



where each scenario reverts to the preceding one if the 
condition on ( or M is not true (e.g., qj = 0.2 when 
M < 2%A/Edd in all three scenarios). In the optimistic 
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Figure 2. Predicted 5 GHz curves for the tidal disruption 
candidate D3-13 (Gezari et al. 2009) for the internal jet 
model in the optimistic scenario (Eq. 3a). The existing up- 
per limits on the 1.4 GHz flux density (Bower 2011) are 
shown with blue triangles, pointing left and right to indicate 
the uncertainty on the time of disruption. The downward 
pointing triangles labeled "2012" show our upper limits to 
the 5 GHz flux. For cos i < f3j, the peak of the light curve 
decreases in time and magnitude as i increases, because a 
fixed observed frequency corresponds to a smaller radius 
where the jet becomes optically thin to synchrotron self- 
absorption (i.e., z ssa cx 5/v). For cosi > f3j the light curve 
is compressed due to time retardation. 



scenario (a), the TD jet behaves like a radio- loud quasar 
at all times. In the most conservative scenario (b), the jet 
becomes radio-loud only when the accretion drops below 
< 2%MEdd (Maccarone 2003). While in c, the system starts 
with a radio-loud burst during the onset of the accretion. 
We consider c the most realistic scenario, since it resembles 
most closely what is observed in X-ray binaries. 

Besides qj (Eq. 3) the internal model requires a jet 
Lorentz factor and the disk luminosity as a function of time. 
The latter is obtained from the fallback rate of the stellar 
debris for a pericenter passage at the disruption radius, 
capped at the Eddington limit. We vary Tj between 5, the 
default value in van Velzen et al. (2011), and Tj — 2. In 
Fig. 2 we show an example light curve of the TDF candi- 
date D3-13 (Gezari et al. 2009) for Tj — 5, using the black 
hole mass as estimated from the luminosity of the host. 

In the internal TD jet model, the typical time scale of 
the light curve is set by the radius where the jet becomes op- 

_l r 2/3 



tically thin to synchrotron self-absorption, 



cx 5v L 



This is different from our off-axis version of the external 
model, where we assumed that all emission happens at the 
head of the jet, which explains the dissimilar scaling of the 
peak of the light curve with the Doppler factor for the in- 
ternal and external model. 



3. Observations 

In Table 1 we summarize the published radio follow-up ob- 
servations of TDF candidates that were discovered at op- 
tical to soft X-ray wavelengths. To increase this sample, 
we selected all TDF candidates with an estimated time of 
disruption after 2004 for follow-up observations. This limit 
is used since the internal model of TD jets is no longer 
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name 


to 


F v 


V 


At 




(yr) 


(mJy) 


(GHz) 


(yr) 


NGC 5905 1 


1999.2 


< 0.15 


8.5 


6.0 


D3-13 2 


2004.5 


< 0.15 


1.4 


1.8 


TDE2 3 


2007.8 


< 0.10 


8.4 


1.1 


CSSf002f7 4 


2010.2 


0.50 


7.9 


0.4 


SDSS Jf20f+30 5 


2010.4 


< 0.22 


4.8 


1.4 



Table 1. Published radio follow-up observations, prior to 
this work, of TDF candidates that were discovered at op- 
tical to soft X-ray frequencies. Except for CSS100217, we 
show the 5-ct upper limits on the radio flux {F v ). At denotes 
the time with respect to the estimated time of disruption 
(to). We note that most of the TDF candidates that were 
discovered in the nineties have post disruption radio upper 
limits from large radio surveys (e.g., NVSS, Condon et al. 
1998), these are listed in Komossa (2002). References: (1: 
Komossa 2002), (2: Gezari et al. 2008; Bower 2011), (3: van 
Velzen et al. 2011), (4: Drake et al. 2011), (5: Saxton et al. 
2012). 



name 


redshift 


M BH 


iint 


<t{F v ) 


At 






Mq x 10 7 


(min) 


(^Jy) 


(yr) 


Dl-9 1 


0.326 


5 


30 


9 


8.0 


D3-13 1 


0.370 


2 


18 


8 


7.6 


TDE1 2 


0.136 


1 


28 


10 


5.4 


D23H-1 3 


0.186 


5 


28 


8 


4.8 


TDE2 2 


0.252 


5 


25 


12 


4.3 


PTFlOiya 4 


0.224 


1 


18 


8 


1.6 


PSl-10jh 5 


0.170 


0.4 


39 


15 


0.71 



Table 2. Jansky VLA observations at 5 GHz of TDF can- 
didates discovered at UV to optical frequencies. We list the 
redshift and estimated black hole mass of these candidates 
in the second and third column, respectively. No significant 
emission was detected at the phase center of the images. 
We list integration time after removal of interference, the 
rms of the images, and the time of the observations with 
respect to the estimated time of disruption. References: (1: 
Gezari et al. 2008), (2: van Velzen et al. 2011), (3: Gezari 
et al. 2009), (4: Cenko et al. 2012a), (5: Gezari et al. 2012). 



valid when the jet slows down significantly. Our sample 
also includes TDFs with existing radio upper limits, since 
the radio emission can be delayed with respect to the time 
of disruption, see previous section. We removed CSS100217 
from the sample because it is detected at 1 GHz before 
the time of disruption with a flat spectral index, indicating 
an AGN origin for the radio emission (Drake et al. 2011). 
SDSS J1201+30 was not selected for follow-up observations 
because this TDF was published after our observations were 
scheduled. Details of the data reduction of the remaining 
six candidates are summarized below. 

The radio observations were carried out on the Karl G. 
Jansky Very Large Array on 29 January 2012 under pro- 
gram 12A-005. We observed at a central frequency of 5.0 
GHz with 16 subbands each with 64 2 MHz channels, span- 
ning 2 GHz of total bandwidth. The VLA was in the C 
configuration yielding typical angular resolution of 4 arc- 
sec. The total observing time was 2.5 hrs, with integration 
times for individual TDF sources varied from 18-30 min. 
Phase calibration was carried out by making short observa- 
tions of nearby point source calibrators every 10 minutes, 



while amplitude and bandpass calibration was achieved us- 
ing an observation of 3C 286 or 3C 48 at the beginning or 
end of each observing run. The data were reduced follow- 
ing standard practice in the Astronomical Image Processing 
System (AIPS) software package. 

In addition to these data, we identified one public data 
set from the VLA archive (project AS 1020) for the TDF 
candidate PSl-lOjh (Gezari et al. 2012). These observation 
were made with the VLA on 29 March 2011 in the B config- 
uration with two subbands (each with 64 2 MHz channels) 
centered at 4.83 and 4.96 GHz, for a total bandwidth of 
256 MHz. The calibration and imaging of these data was 
similar to the method described above. 

Our final sample that we shall use to constrain TD jet 
models thus consists of seven TDF candidates that were 
observed with the Jansky VLA. We summarize the results 
of these observations in Table 2. 



4. Analysis 

In this section we first compute the constraints that can be 
placed on TD jet models using our Jansky VLA follow-up 
observation and then consider the potential of radio tran- 
sient surveys. 

4.1. Constraints from follow-up observation 

If we assume that the angle between the observer and the 
jet is drawn from a uniform distribution (on a sphere), we 
can calculate the probability of non-detections for a given 
flux density limit. One simply has to find the largest angle 
for which the predicted flux is above the flux limit and then 
calculate the probability to observe a jet within this angle. 
The flux limit is set at twice the rms of the radio image 
of each TDF. (This is lower than the limit for a blind- 
detection experiment since we use the threshold to find the 
probability of a non-detection, not to claim a discovery.) In 
Table 3 we list the results of this exercise. 

The probability that all seven TDFs in our sample 
hosted jets, but were not detected due to Doppler boosting 
is P-j — HiPi, with Pi being the probability of the observa- 
tions of each TDF candidate, as listed in Table 3. We also 
consider the possibility that, given our observations, at least 
one of the seven TDFs hosted a jet, P>\. This is obtained by 
taking the mean value of the product of all combinations 
of the seven Pj's (e.g., the probability that only one jet 
was launched is P\ = JV Pi/7)- For the optimistic scenario 
of the internal model (sec. 2.2) with Tj = 5, four of the 
seven TDF candidates (D23H-1, TDE2, PTFlOiya, PSl- 
lOjh) should have yielded a detection above the 2-a level, 
hence Pf = 0, while P>i=2%. The probability that all of 
the other three TDF candidates hosted a jet is 1.7%. For 
the most conservative scenario (Eq. 36), P-j — 48%, while 
for the realistic scenario (Eq. 3c) this is lower at 21%. In 
Fig. 3 we show P-j and P>\ for lower Lorentz factors; at 
Tj < 3, the hypothesis that all seven TDF hosted a jet is 
ruled out at 95% confidence for all three scenarios of the 
internal jet model. 

Our upper limits also constrain the possibility that a jet 
similar to Sw 1644+57 was launched after the disruption. 
To place the Jansky VLA observations on the estimated off- 
axis light curve (Eq. 2), we equate the time of disruption 
to the time of the Swift trigger and we scale the flux using 
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name Internal jet model Sw f 644+57, off-axis 





a 


c 


b 


r 3 = 2 


r, oc t~ - 2 


Df-9 


39 


78 


83 


49 


46 


D3-13 


62 


89 


91 


52 


50 


TDEf 


7 


92 


100 








D23H-f 





52 


70 








TDE2 





75 


98 


20 


19 


PTFlOiya 





86 


95 








PSf-fOjh 





95 


97 









Table 3. Probability (%) that the jet orientation is such 
that the predicted flux is below the 2cr-level of the 5 GHz 
observation. Zero probability implies the predicted flux is 
above the threshold even for i 2 = 7r/2, while Pi = 100% 
implies the data cannot constrain the model. In the sec- 
ond to third column we list the results for the internal jet 
model, for the optimistic to the conservative scenario (Eq. 
3), using Tj = 5. In Fig. 3 we show the results for lower 
Lorentz factors. In the fourth and fifth column we give the 
probably of detecting a jet that is identical to Sw 1644+57, 
but observed off-axis, using two different estimates of the 
light curve past the last available observation (see sec. 2.1). 



3 0.3 



o.o 



Always radio-loud (a) 

Conservative (b) 

^™ Realistic (c) 


•>* 


✓ 







2.0 2.5 3.0 3.5 4.0 4.5 5.0 



r, 

Figure 3. The probability of our data (i.e., no flux above 
the two times the image rms) for the three scenarios of 
the internal model (Eq. 3). The solid lines show P7 the 
probability that all seven TDFs we observed indeed hosted 
a jet; the dashed lines show P>i the probability that at 
least one hosted a jet. For Tj < 3 the former hypothesis is 
ruled at 95% confidence for all scenarios. 



(dL,Sw/d,L) 2 , with d,L,s-w being the luminosity distance of 
Sw 1644+57. From Fig. 1 we see that our upper limits on 
the radio flux of four TDFs (TDE1, D23H-1, PTFlOiya, and 
PSl-lOjh) are inconsistent with the estimated off-axis light 
curve of Sw 1644+57 for all viewing angles. The probability 
that the other three flares hosted a radio transient identical 
to Sw 1644+57 is about 5% (for both versions of the late- 
time evolution we considered in sec. 2.1). 

4.2. Constraints from (future) radio transients surveys 

A different method to test whether jets like Sw 1644+57 are 
common to stellar tidal disruptions is to compute the rate 
of these transients. The snapshot rate (or areal density) at 
a given flux density limit -Fl/.lim can be estimated directly 
from the Sw 1644+57 light curve: 

i?(i^i im )~8xio- 3 iy 2 (§^y /2 

AT ^ TD d 2 (4) 

10- 5 5 x IO- 3 Mpc- 3 W 

Here AT is the time in years that the light curve of 
Sw 1644+57 is above -F„Sw) Ntd is the stellar tidal disrup- 
tion rate, and pbh is the black hole density. If Sw 1644+57 
was a typical stellar tidal disruption, this rate should be of 
the same order as the TDF rate inferred from soft X-rays 
(Donley ct al. 2002) or optical (van Velzen and Farrar 2012, 
in prep) surveys, i.e., Ntb ~ 10 yr. 

The 5 GHz light curve of Sw 1644+57 implies AT w 1 yr 
for T^sw = 20 mJy. For Tj — 2, we can thus obtain the 
snapshot rate for a radio variability survey with a threshold 
at 10 mJy: i?(10mJy) = 5x 10~ 3 deg -2 . This rate is close 
to the existing upper limits on the snapshot rate at 5 GHz 
(e.g., Scott 1996; Bower et al. 2007, 2011) - see Frail et al. 
(2012) for a review. So near-future radio variability survey 
will either measure or constrain A/td for transients similar 
to Sw 1644+57. 



5. Conclusion & Discussion 

We obtained upper limits at the ~ 10 \i Jy level of the 5 GHz 
flux of seven stellar tidal disruptions events that were dis- 
covered with optical/UV imaging surveys. This is three or- 
ders of magnitude lower than the recently discovered TDFs 
with radio emission. Suggesting that stellar tidal disrup- 
tions come in different flavors, ranging from radio-loud to 
radio-quiet (or radio-silent). To explore how this conclusion 
would be biased by the large possibile parameter range in- 
herent to TDFs, we compared our upper limits to currently 
available jet models, taking into account Doppler boosting 
and temporal evolution of the radio emission. 

We used our observations to constrain the jet model of 
van Velzen et al. (2011). For a jet Lorentz factor of Tj = 5, 
we can rule out the optimistic ( "alway radio-loud" ) scenario 
for four of the seven flares. The probability that the other 
three TDF candidates did launch such jets, but are not 
detected because Doppler boosting reduced the flux below 
two times the image rms is only 4%. The hypothesis that all 
events hosted jets that only becomes radio-loud when the 
fallback rate drops below 2% of the Eddington accretion 
rate (i.e., as observed in stellar mass black holes) is less 
constrained. Only for jets with Tj < 3 this hypothesis is 
ruled out at 95% CL. 

We have also investigated the possibility that our sam- 
ple of TDFs hosted a jet which is identical to Sw 1644+57, 
but oriented at a larger angle between the observer and 
the jet. For four of the seven flares, the estimated off-axis 
light curves of this relativistic TDF are inconsistent with 
the non-detection for all possible observer angles. The hy- 
pothesis that all of the other TDFs hosted jets identical to 
Sw 1644+57 is ruled out at the 95% confidence level. These 
results are not sensitive to our assumption that the time of 
disruption equals the time of the Swift trigger. If the hard 
X-rays of the jet are emitted only after ten times the fall- 
back time (~ 1 yr), the predicted off-axis flux is increased 
by just 50%. A more serious caveat is that the off-axis light 
curves we used in this work are only valid for circumnuclear 
environments that are identical to the host of Sw 1644+57. 
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This is not likely to be the case: the blue colors of opti- 
cal/UV flares imply little optical extinction (i.e., redden- 
ing), while for Sw 1644+57 this extinction is much higher, 
Ay = 3 — 5 mag (Bloom et al. 2011). Finally we note that 
the black hole mass of Sw 1644+57 may be a factor 5-10 
smaller than the median black holes mass of thermal TDFs 
that we followed-up. For a pericenter passage at the disrup- 
tion radius, this implies that for the flares in our sample, the 
duration of the super-Eddington fallback rate is expected 
to be about a year shorter and the total jet energy could 
be factor 3 higher (for an accretion rate that is capped at 
the Eddington limit). A more sophisticated treatment of 
the off-axis light curves in the external model should take 
these differences into account, e.g., using radiative trans- 
fer onto the output of 2D hydrodynamical simulations (van 
Eerten & Wijers 2009) for a range of black hole masses and 
environments. 

Definite proof that relativistic TDFs with evidence for 
jetted emission are an intrinsically different class can be 
obtained by radio transient surveys. For a disruption rate 
that is of the same of order as the rate of thermal TDFs, the 
areal density of radio transients like Sw 1644+57 (Eq. 4) 
almost exceeds the current upper limits. Near-future radio 
variability surveys, such as VAST (Murphy et al. 2012), 
ThunderKAT, which is part of or MeerKAT (Booth et al. 
2009), or LOFAR 2 , will either detect tens to hundreds of 
TD jets per year, or conclude that the rate of stellar tidal 
disruptions with jets is lower than the rate of (radio-quiet) 
thermal TDFs. 

If a division between tidal disruptions with and without 
jets indeed exists, it presents a challenge to the idea that 
radio-loudness can be explained by state changes of accre- 
tion disk. Some authors have argued that the spin of the 
black hole is important to produce stellar tidal disruption 
jets (e.g., Lei & Zhang 2011; Krolik & Piran 2011; Stone & 
Loeb 2012). Testing this idea requires observations of the 
disk emission of tidal disruptions with jets to show that 
this component is similar to TDFs without jets. Otherwise, 
differences in the evolution of the stellar debris, which may 
depend on the orbit of the star or the presence of a pre- 
existing accretion disk, can also explain the observed range 
in radio-loudness. 
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